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O 

' We study the sign problem in lattice field theory with a 6 term. We apply the maximum 

I entropy method (MEM) to flattening phenomenon of the free energy density f{9), which 

K , originates from the sign problem. In our previous paper, we applied the MEM by employing 

^ ' the Gaussian topological charge distribution P{Q) as mock data. In the present paper, we 

consider models in which 'true' flattening of f{6) occurs. These may be regarded as good 
examples for studying whether the MEM could correctly detect non trivial phase structure. 
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§1. Introduction 

A 9 term in QCD is expected to cause interesting dynamics. ^^'^^ Lattice field 
' theory may be a suitable tool as a non perturbative method to study its effect. ^^"''^ 

However, it suffers from the complex action problem or the sign problem, when a 6 
{J^ \ term is included. A conventional technique to avoid the problem is to perform the 

<^ ' Fourier transform of the topological charge distribution P{Q), which is calculated 

with real positive Boltzmann weight, and to calculate the partition function Z(6).^^'^^ 
However, this still causes flattening of the free energy density f{9) originated from 
p) ' the error in P{Q)}^^'^^^ This phenomenon misleads into a fictitious signal of a first 

1^ . order phase transition. Some approaches have been tried to overcome the sign 

' problem. ^^^'^'^^ In our previous paper^^^ (referred to as (I) hereafter), we applied to 

this issue the maximum entropy method (MEM), which is suitable to deal with so- 
^ . called ill-posed problems, where the number of parameters to be determined is much 

^ I larger than the number of data points. The MEM is based upon Bayes' theorem. It 

derives the most probable parameters by utilizing data sets and our knowledge about 
these parameters in terms of the probability. The probability distribution, 
which is called posterior probability, is given by the product of the likelihood function 
and the prior probability. The latter is represented by the Shannon-Jaynes entropy, 
which plays an important role to guarantee the uniqueness of the solution, and the 
former is given by x^- The task is to explore the image such that the posterior 
probability is maximized. 

In the paper (I), we used the Gaussian P{Q) and tested whether the MEM 
would be effective, since the Gaussian P{Q) is analytically Fourier-transformed to 
Z{9). For the analysis, we used mock data by adding Gaussian noise to P{Q). The 
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results we obtained are the followings; In the case without flattening, the results of 
the MEM agree with those of the Fourier transform method and thus reproduce the 
exact results. In the case with flattening, the MEM gives smoother f{9) than that 
of the Fourier transform. Among various default models investigated, some images 
with the least errors do not show flattening. 

If the reason for the disappearance of flattening in the above analysis is the 
smoothing of the data, which is the characteristic of the parameter inference such 
as the maximal likelihood method, one may wonder if the MEM may not be a 
suitable method for detecting singular behaviors such as a phase transition. In order 
to investigate this issue, we study whether the MEM would be applicable to 'true' 
flattening of f{9) or singular behaviors like a first order phase transition occurring 
at a finite value of it). By 'true', it is meant that fiattening is not originated 
from the error of P{Q) but occurs due to the data themselves. For this purpose, we 
consider a simple model based on mock data. It is a model which consists of the 
Gaussian P{Q) except at Q = 0. For Q = 0, a constant is added to the Gaussian 
P{Q). This additional contribution at Q = causes the fiattening behavior in the 
large 9 region (6* < tt). We apply the MEM to this model and study whether true 
flattening would be properly reproduced. 

Although this model develops true flattening, it does not reveal the singularity in 
f{9) at flnite value of 9[^ vr). In order to study such a singular behavior mimicking 
a first order phase transition, we consider a model which utilizes P{Q) obtained by 
the inverse Fourier transform of a singular f{9). In order for the singularity and 
flattening to be compatible in such a model, obtained P{Q) oscillates and could take 
negative values for large values of Q. Although P{Q) looses the physical meaning 
as the topological charge distribution for their negative values, we regard this as a 
mathematical model to study the singular behavior in terms of the MEM. 

We apply the MEM to these models and study (i) whether true flattening would 
be properly reproduced and (ii) how the errors in P{Q) are related with the results 
of the MEM and the numerical Fourier transform method. Although these models 
are simple, we expect that this analysis would serve to obtain knowledge if the MEM 
could be utilized to study the phase structure, when it were applied to other models 
suffering from the sign problem, say, flnite density QCD. 

This paper is organized as follows. In the following section, we give an overview of 
the origin of flattening and a sketch of the MEM. A model developing true flattening 
is also introduced. In § |31 the MEM is applied to this model. In § ^ another model 
developing a singular behavior is introduced and analyzed. A summary is given in 

§2. MEM and a model 

2.1. Flattening of f {9) 

Let us briefly explain flattening of f{9). The partition function Z{9) can be 
obtained by Fourier-transforming the topological charge distribution P{Q), which is 
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calculated with real positive Boltzmann weight. 



/ [dzdz] 



-S+i9Q{z,z) 



where 

The measure [dzdzjq in Eq. ()2-2p represents that the integral is restricted to config- 
urations of the field z with the topological charge Q, and S denotes an action. 

The distribution P{Q) obtained from Monte Carlo simulations can be decom- 
posed into two parts, a true value P{Q) and an error of P{Q). When the error at 
(5 = dominates, the free energy density becomes 

f{e) = -^\ogZ{e) « -llog{ j;p(Q)e^^Q + 5P(0)} 



1 

v 



log{e-^^» + 5P(0)}, (2-3) 



where f{6) is the true free energy density and (^P(O) is the error at Q = 0. The 
value of e"^-^*-^) decreases rapidly as the volume V and/or 6 increases. When the 
magnitude of e"^-^*^^) becomes comparable to that of (5P(0) at 9 9f, f{6) becomes 
almost constant, f{9) ~ const, for 6 > 9f. This abrupt change of the behavior is 
misleadingly identified as a first order phase transition at 6 k, 9{. This is (fictitious) 
flattening phenomenon of f{6). The error in P{Q) could also lead to negative values 
of Z{9). We refer to this, too, as flattening, because its origin is the same. In the 
case of a large volume, an exponentially increasing amount of data is needed. 

2.2. MEM 

In this subsection, we briefly explain the MEM in the context of the 9 term. See 
(I) for details. Instead of dealing with Eq. (|2-ip . we consider 

/TT p-ieq 
J9^^Z{9). (2-4) 

The MEM is based on Bayes' theorem. In order to calculate Z{9), we maximize the 
probability 

prob(Z(«)|F(Q)./) =prob(ZW|/)E;;^i^l|m^ (2.5) 

where prob(j4) is the probability that an event A occurs, and proh{A\B) is the 
conditional probability that A occurs under the condition that B occurs. 
The likelihood function proh{P{Q)\Z{9), I) is given by 

-1^2 

prob(P(Q)|Z(^),/) = ^-^, (2-6) 
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where Xl is a normalization constant, and is represented by 

^ - PmCg^^'iP'-^HQ') - PiQ'))- (2-7) 

Q,Q' 

Here, P^^\Q) is constructed from Z{9) through Eq. l\2-4\\ . and P{Q) denotes the 
average of a data set {P{Q)}; 

1 ^'^ 

P{Q) = j^Y.P^'^(Q^' (2-8) 

1=1 

where A'^^ represents the number of sets of data. The matrix represents the 
inverse covariance matrix obtained by the data set {P(Q)}. 

The prior probabihty pioh{Z{6)\I) is represented in terms of the entropy S as 

prob(Z(^)|/,a) = -^, (2-9) 
Xs[a) 

where a is a real positive parameter and Xs{c() denotes an a-dependent normaliza- 
tion constant. With regard to S, the Shannon-Jaynes entropy is employed: 



S= dO 



z{e) - m{e) - z{e)\og ^^^^ 



(2-10) 



m{e) 

where m{9) is called the 'default model'. The posterior probability thus amounts to 

prob(Z(0)|P(Q), /, a, m) = ^ ^ ^ V (^-ll) 

where it is explicitly shown that a and m are regarded as new prior knowledge in 
pvoh{Z{e)\P{Q),I,a,m). 

For the information /, we impose the criterion 

Z{e) > (2-12) 

so that pvoh{Z{e) < 0\I,m) = 0. 

We employ the following procedure for the analysis. Details are given in 

(I). 

1. Maximizing W to obtain Z^^'^O) for a fixed a: 

In order to find the image in functional space of Z{9) for a given a, we 
maximize 

W^-^X^ + aS (2-13) 



so that 



^ i-y + c^S) =0. (2-14) 



5Z{ey 2 

The parameter a plays the role determining the relative weights of S. 
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In the numerical analysis, the continuous function Z{9) is discretized; 

Z{9) Z{9n) = Zn- The integral over 9 in Eq. (|2-4|1 is converted to the finite 

summation over 9; 



^ 27r 

n=l 

E 

n=l 



COS tin] 

vr 



(j = 1) 



Zn (otherwise) 



(2-15) 



?i=i 



where j = 1, 2, • • • ,Nq and n = 1, 2, • • • , Ng. Note that Pj denotes P{Q) at 
Q = j — 1 and A'^g < A'^g. Here we used the fact that P{Q) and Z{9) are even 
functions of Q and 0, respectively. 

2. Averaging -E^"^: 

The Qf-independent final image, denoted as Zn or Z{9), can be calculated 

r(«) 



by averaging the image Zn according to the probability. 

Zn = j da prob(a|P(Q),I,m)Z^"). 
The probability prob(Q|P(Q), /, m) is given by 
prob(a|P((5), /, m) = — — exp 

where is a normalization constant, and A{a) = ^Ylk a+Afc • Here the 
values Afc's are eigenvalues of the real symmetric matrix in 9 space; 



(2-16) 



(2-17) 



dZfndZn 



Zn 



(2-18) 



In averaging over a, we determine a range of a so that pioh{a\P{Q), I,m) > 
^ xprob(a|P((5), "^) holds, where pToh{a\P{Q), I, m) is maximized at a = d. 



3. Error estimation:-'^^)'^'^) 

The error of the final output image Zn is calculated based on the uncertainty 
of the image, which takes into account the correlations of the images Zn among 
different values of n.^^) 



2.3. A model 

In (I), we used the Gaussian P{Q) for the MEM analysis. We parameterized the 
Gaussian P{Q) as 

PG{Q)^AeM-^Q\ (2-19) 

where, e.g., in the case of the 2-d U(l) gauge model, c is a constant depending on 
the inverse coupling constant /?, and V is the lattice volume. Here V is regarded as a 
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parameter and varied in the analysis. The constant A is fixed so that YIq PciQ) = 
1. In order to study the question whether the MEM could be a useful tool for 
reproducing true flattening, we consider a simple model. 

Suppose in some lattice theory that calculated P{Q) were slightly deviated from 
the Gaussian one only at Q = 0: 

P{Q) = PG{Q)+Poy~5Qo^ (2-20) 

where pq is a positive constant and 6qq is the Kronecker symbol. The partition 
function is analytically obtained by use of the Poisson sum formula as 

+ Po}, (2-21) 

where A' is the normalization constant, which is determined so that -2^pois(0) = 1- 
Because the first term in Eq. ()2-21|) is monotonically decreasing function of ^, the 
role of the constant term is to cause a flat distribution at large values of ^ < vr). 
Let us call this model a flattening model. We use the distributions P{Q), Eq. (|2-2U|) . 
to generate mock data by adding Gaussian noises to them as described below. 

2.4. Mock data 

For preparing the mock data, we add noise with the variance 5x P{Q) to P{Q). 
This is based on the observations in the procedure which was employed to calculate 
P{Q) in the simulations of models such as the CP^~^ model. This turns out to 
yield a ratio of the error to the data to be almost constant; \5P{Q) / P{Q)\ ~ const. 
The parameter 6 directly affects the covariance matrix, and indirectly influences the 
error 6Z{6) of the image, which is the uncertainty of the image. As done in (I), 5 is 
mainly fixed to 1/400. In order to investigate the effects of 6 on the final image, we 
also use its smaller values. 

A set of data consists of P{Q) with the errors for Q = to A^g — 1. Employing 
such sets of data, we calculate the covariance matrices in Ea. (|2-7() with the jackknife 
method as 

-, A^d _ 

Cq,Q' = ^ (pffl(Q)-P(Q))(p(')(QO-mO). (2-22) 

where P^^\Q) denotes the /-th data of the topological charge distribution and P{Q) 
is the average Ea. ()2-8|1 . The inverse covariance matrix is calculated with precision 
that the product of the covariance matrix and its inverse has off-diagonal elements 
at most 0(10-23). We take Nd = 30 as in (I). 

The number of degrees of freedom in 6 space, Ng, is taken to be larger than 
that of the topological charge Ng: Nq is set to be 28.^^^ The number Ng is chosen 
so as to satisfy P{Q) > 10^^'^ and to satisfy the quality of the covariance matrix 
described above. It varies depending on which type of model is used and on V. Note 
that in order to reproduce 2(9) which ranges many orders, the analysis must be 
performed with quadruple precision. In order to find a unique solution in the huge 
configuration space of the image, we employ the singular value decomposition. 



^pois(^) = A' 



InV 
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As to the default model m{9) in Eq. (|2-in|) . we study various cases: (i) m{9) = 
const, and (ii) m{9) = exp(— ^24^7^^). In the case of (i), we show the results of 
m{6) = 1.0. The case (ii) is the Gaussian default, where the parameter 7 is varied 
in the analysis. 

In the next section, we investigate the following questions: (i) Does the MEM 
reproduce true flattening? (ii) How does the final image of the partition function 
depend on pQ as well as on the magnitude of the error in P{Q)'! 



3.1. Flattening model 



§3. Results 
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Fig. 1. Partition function Z{9) calculated by numerical Fourier transform for V = 50. Data for 
P{Q) in Eq. ll^^ are used with po = 4.5 x 10-^ 4.5 x 10'^, 2.0 x 10"^ and 0. Here, 5 = 1/400. 
Clear flattening is observed in each case. The result for po = looks similar to that for 2.0 x 10~'^ . 
The latter indicates somewhat smaller errors than the former. 

When noise is added, the partition function for Eq. (j2-2()j) in the flattening model 
is given by 

{EQ^G(Q)e^^« +P0 + EQ^Qy"^} 



m 



B 



(3-1) 



where B = '}2q{Pg{Q) +Po'5qo + ^P{Q)^- Figure ^ displays the behavior of Z{9) 
calculated by the numerical Fourier transform of Eq. (|2-2Up for = 50 and for 
various values of poi Po = 4.5 x 10~^, 4.5 x 10^'^ and 2.0 x 10~'^. The parameter c 
in Eq. (|2-19|) is fixed to 7.42 throughout the paper. The parameter 6 is chosen to be 
1/400. FigiureQalso displays Z{6) for po = 0, which was studied in (I) in detail. 
Clear flattening is observed in each case. Note that the behavior of Z{6) and its 
error are similar for the po = 2.0 x 10^^ and = cases. In order to understand 
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Fig. 2. Z^"\9) obtained in the procedure 1, Eq. II2-14I I. for various values of a in the flattening 
model. Here, V = 50, po = 4.5 x 10~^ and S = 1/400. The default model is chosen to be the 
constant m{9) — 1.0. 
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Fig. 3. Probability prob(a|P((5), /, m) = P(a) for V = 50, po = 4.5 x 10"^ and S = 1/400 in the 
flattening model. The default model is chosen to be the constant m{9) — 1.0. 



what the origin of these flattening is, let us investigate Eq. (|3-1() briefly: When pQ is 
so chosen that for small values of 9, the first term dominates in the numerator, while 
for large values of 9, the constant term dominates, the resultant 2(9) behaves like 
the partition function with flattening. For the latter case, Eq. is approximated 
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Fig. 4. .2(0) for F = 50, po = 4.5 x 10~^ and 5 = 1/400 in tiie flattening model. The default model 
is chosen to be the constant m{9) — 1.0. Error bars are also displayed but invisible. 



by a constant 

Z{9) « |. (3-2) 
For V = 50 and po = 4.5 x 10^^ , B is numerically estimated as 4.45 and thus 

m ^ = 1-0 X 10-2. (3.3) 
4.45 

For other choice of pq, 4.5 x 10"^ and 2.0 x lO^"^, po/B becomes 1.0 x 10^'^ and 
4.5 X 10~^, respectively. These values agree with those of flattening observed in 
Fig. n and we thus find that these flattening are not caused by the error in P{Q) 
but are due to the additional term in Eq. (|2-20() . On the other hand, to flattening 
for the Po = case in Fig. Eq. (|3-2|) does not apply obviously. What is observed is 
not caused by the constant term but comes from the error in P{Q) as described in 
section [TTl The former type of flattening is true flattening, because its origin lies in 
data of P{Q), while the latter {pq = 0) is fictitious one, because it comes from the 
error in P{Q). 

Having distinguished two kinds of flattening, we investigate whether the MEM 
can properly reproduce the true one. We employ the case for pQ = 4.5 x 10"^, as 
an example. Figure HI displays Z^^^O), which is calculated by Eq. H2-14() . for various 
values of a. The default model is chosen to be the constant, m{6) = 1.0. We 
find that 2^°'^ (9) are independent of the values of a. They show clear flattening in 
agreement with the result of the Fourier transform and also with that of the Poisson 
sum formula, Eq. (|2-21|) . 

In order to calculate the final image of the partition function, we need the proba- 
bility prob(a|P((5), /, m) as shown in Eq. H2-16() . For m{9) = 1.0, prob(a|P((5), /, m) 
is plotted in Fig. |31 A peak of prob(a|P((5), /, m) is located at a ~ 12. The inte- 
gration to obtain the final image Z{0) thus becomes trivially simple because of the 
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Fig. 5. Z{e) for 1/ = 50, po = 1-0 x 10"^ and 5 = 1/400 in tiie flattening model. The default 
model is chosen to be the constant m{9) — 1.0. Error bars are also displayed but invisible for 
the MEM results. The Fourier method yields negative values for Z{6) a.t 6 > 2.5. 



Table L The final image Z{0) in the flattening model at = 2.32 and 9 — 3.07 for various values 
of Po- The default model is chosen to be m{0) = 1.0. The exact values 2pois(6') are also listed. 
The results of the Fourier method are denoted as Zfom- 



Po 


^pois (2.32) 


-Efour(2.32) 


^(2.32) 


4.5x10"^ 
4.5x10"^ 
1.0x10"^ 


9.801x10"^ 
1.094x10"^ 
3.342x10"" 


9.6(l)xl0-^ 
9.1(1.5)xl0-" 
1.5(1.5) X 10"" 


9.66(l)xl0~^ 
8.99(2)xl0-" 
1.78(8)xl0"^ 


Po 


2pois(3.07) 


2four(3.07) 


2(3.07) 


4.5x10"^ 
4.5x10"^ 
1.0x10"^ 


9.686x10"' 
9.772x10"" 
2.174x10"" 


9.5(2) X 10"^ 
7.5(1.7)xl0"" 
-1.5(18.0)xl0"5 


9.521(2)xl0"^ 
7.410(6) X 10"" 
3.6(2) xlO"« 



fact that Z^'^\9) is independent of a in the relevant range for the integration of the 
procedure 2. Thus Z{9) agrees with Z^'^\9) in Fig. |2J and the MEM yields true 
flattening. Note that this situation is unchanged when the Gaussian default model 
is used, where 7 is varied between 0.1 and 3.0. 

With other choice of the value of po, we repeat the same procedure. For po = 
4.5 X 10~^, the MEM yields true flattening for the final image Z{9) as plotted in 
Fig. 01 The errors are also included but invisible in the figure. The result of the MEM 
agrees with the behavior of the exact partition function (Poisson). For smaller values 
of pq, however, flattening is not reproduced correctly. This is displayed in Fig. [HI for 
Po = 1.0 X 10"'^; although Z(9) in the MEM flattens for large values of 9, its height is 
about 4 times smaller than the exact one (Poisson). This analysis shows that there 
is a critical value of po above which the MEM is effective, and it is found that pQ 
~ 2.0 X 10"^ in the present case. This critical value is associated with the magnitude 
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of the error in P{Q) or 6. This will be discussed later. The result of the Fourier 
method, Zfour(^)> is also shown as a comparison in Fig. |S1 Here, ZiomiO) receives 
large errors for 2.3 ^ < 2.5. For 6 > 2.5, Zfo^^^{9) takes negative values and are not 
plotted in the figure. Table Q] displays the numerical values of Z(9) at two typical 
values of 6, 2.32 as an intermediate region and 3.07 as the one near tt, for three 
values of po- The exact values -Zpois(^) and Ziouvid) ai'e also listed. 

3.2. 6 dependence 

In the model discussed above, the applicability of the MEM is deeply associated 
with the magnitude of the error in P{Q). In order to investigate this aspect, we 
vary the value of 5, which controls the magnitude of the error in the mock data. 
Although it would not be realistic to use data with too small error corresponding to 
huge amount of statistics, it would still be meaningful to investigate the limitations 
of the MEM analysis. This is beneficial of using mock data. 
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Fig. 6. Z{6) in the flattening model obtained by the numerical Fourier transform and the MEM 
for \/ = 50, po = 1-0 X 10~^. The value of 5 is chosen to be 1/400 and 1/600 for a comparison. 
Note that 2four(^) takes negative values for 6 > 2.5 in the case S = 1/400 and thus these data 
are not shown. Error bars for Zfourid) with 5 = 1/600 are large, while those of the MEM are 
invisible. 

For the flattening model, we fix the value of po and vary that of 6. Let us take 
an example of po = 1-0 x 10~^, which is about half the critical value of pq for the 
6 = 1/400 case. As already seen in Fig. both of the methods were not effective 
in reproducing flattening. We vary the value of 6 from 1/400 down to 1.0 x 10~^. 
It is found that already for 5 = 1/500, the MEM reproduces flattening reasonably 
well. Figure El displays the results of the MEM analysis and those of the numerical 
Fourier transform for 5 = 1/400 and 1/600 as a comparison. The error bars shown 
in the figure are those of Fourier transform, while those for the MEM are invisible. 
For 5 = 1/600, although the errors for the Fourier method are large at ~ vr, the 
central values of Zfouri^^) are approximately equal to those of the MEM, in contrast 
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to the 6 = 1/400 case where ^four(^) takes negative values for 9 > 2.5. At 6 = 3.07, 
the relative error of Z{our{(^) goes down to 50 % for 5 = 1.0 x and 4.6 % for 
<5 = 1.0 X 10"^ while it is about 100% for 6 = 1/600. 

When po is chosen to be smaller than 1.0 x 10~^, the MEM can recover correct 
images by taking smaller value of 6 than 1/600. For a given value of pq, the MEM 
becomes effective in reproducing correct images, when the value of 6 is reduced to 
the one corresponding to the precision that |5P((5)| ~ po. With this precision 
of data, 2^four(^) still receives much larger errors than those of the MEM. 



§4. Singular f{0) model 



4.1. A mathematical model 

Although the flattening model develops true flattening, it does not reveal a 
singularity in f{9). In order to study such a singular behavior mimicking a first 
order phase transition at finite value of vr), we consider a model which utilizes 
P{Q) obtained by the inverse Fourier transform of a singular f{0). Consider a 
partition function 

^ ' 1 Zi {0i<e<7r), ^^^^ 

where -Zpois(^) is the one obtained from the Poisson sum formula Eq. (|2-21j) with 
Po = 0, and Z{ is a positive constant (= 2pois(^f))- This partition function is 
inverse- Fourier-transformed to obtain P{Q): 



P{Q) = - 

TT 







'f 

deZpoisio)coseQ+ I deZfcosOQ 



(4-2) 



The singular behavior at 9f forces P{Q) to oscillate for large values of Q, and P{Q) 
could take negative values in such a region. Although P{Q) looses the physical 
meaning as a topological charge distribution for their negative values, we regard 
such a set of data of P{Q) as a mathematical model. What we are concerned with 
is rather whether the singular behavior of Z{6) is reproduced. 

In the numerical calculation, however, the singularity appears in a approximated 
way. Also, P{Q) is truncated at an appropriate value of Q satisfying a constraint 
to keep the precision of the covariance matrix, which was described in the subsec- 
tion [231 As a consequence, we consider the case where only positive P{Q) are used 
as well as the one where negative contribution is also taken into account. We use 
the distributions P{Q), Eq. H4-2|) . to generate mock data by adding Gaussian noises 
to them as in the flattening model. We refer to this as singular f{9) model. 

4.2. Results 

We are interested in the behavior at Of and flattening for 9 > Of. We fix Of 
at an appropriate value, say, Of = Tr/2 and vary the value of V. Figure [71 displays 
the result of Fourier transform and the final image Z(0) by the MEM analysis for 
V = 50 and 100. Here, we truncate P{Q) at some value of Q so that P{Q) becomes 
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Fig. 7. Z{6) obtained by tlie numerical Fourier transform and by the MEM for Of = 7r/2 in the 
singular f{6) model. The upper (lower) one corresponds to V = 50(100). The default model is 
chosen to be the constant m{6) = 1.0. Here, 5 = 1/400. Error bars are displayed but invisible 
in both of the methods. Arrows indicate the expected heights of flattening. 
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Fig. 8. Z{6) obtained by the numerical Fourier transform and by the MEM for V — 120 and for 
9f = n/2 in the singular f{6) model. Here, 5 = 1/400. The results with three different default 
models are shown. Error bars are of the Fourier method, while those of the MEM are invisible. 
An arrow indicates the expected height of flattening. 



positive (Q=7 (12) for V = 50(100)). For V = 50, a somewhat acute bend at 0f 
and flattening for 9 ^ Of are observed. Here we used m{9) = 1.0 as a default model. 
Tlie lieight of flattening is in agreement with the input value of Zf = 2pois(7r/2). 
For V = 100, singular behavior at 6{ becomes clearer, and the height of flattening 
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Fig. 9. Z{d) in the singular f{9) model obtained by the numerical Fourier transform and the MEM 
for V — 150 and Of = tt/2. The results for only two values of 5 are plotted: S — 1.0 x 10"'' and 
2.0 X 10~^. The constant default model m(9) = 1.0 is used for the MEM analysis. Error bars 
of Zioui{0) for 5 = 1.0 X 10"'^ are not shown, while those for the MEM are shown but invisible. 



decreases in agreement with its input value. As V increases, the bend at 9i becomes 
more prominent, but the flattening behavior becomes gradually distorted. This is 
because maximal value of Q used for the analysis is restricted in order to keep the 
necessary precision of the covariance matrix as stated in subsection 12.41 ^^^ These 
features can also be seen in the case V = 120 as shown in Fig. |H1 Here, the height of 
flattening agrees with the input value indicated by an arrow in the figure. Figure|Sl 
displays default model dependence of the MEM results, where the default models 
with m{9) = 1.0 and the Gaussian default with 7 = 3.0 and 7 = 5.0 are chosen. For 
two values of 6, the values of 2{0) are listed in Table Ull It is found that the final 
images of Z(9) are independent of these default models, and the magnitudes of the 
errors are almost the same, 1— 3 %. Figure IS] also displays the result of the Fourier 
method Zfouv, where quite large errors are generated for this value of V. In the case 
ofV = 120, P{Q) is truncated at Q = 15. 

As V increases further, the height of flattening decreases, and the image is much 
more affected by the error of P{Q). Here also, the error of P{Q) puts limitations 
to the applicability of the MEM as in the flattening model. For 9^ = it/2, it is 
found that the MEM cannot be applied already aX V = 150, where the height of 
flattening of the image Z[9) comes out about one order larger than the input value 
Zi, although an acute behavior is still seen at Of. This will be discussed in the 
following subsection. 

4.3. 5 dependence 

By taking advantage of the mock data, we can take very small values of 5. When 
the errors in P{Q) are much reduced, the Fourier method provides the results with 
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Table II. Z{6) in the singular f(9) model a.t 6 = 2.32 and 9 — 3.07 for various default models. 
Here, V = 120. 



e 




.27=3.0 


^7=5.0 


^four 


2.32 


1.39(3) X 10-"* 


1.47(3) X 10"'' 


1.46(3) X 10"'' 


1.4 (1.1) X 10"" 


3.07 


5.31(5)xl0"^ 


5.78(5) X 10"^ 


4.65(5) X 10-'^ 


4.7(15.7) X 10"^ 



small errors. Let us take an example oiV = 150 in the singular f{6) model, to which 
both of the methods were not effective for 6 = 1/400 as stated above. We choose 
6 = 2.0 X 10~^ and show the results in Fig. 1^1 Those for 6 = 1.0 x 10~^ are also 
plotted as a comparison. It is found for 6 = 1.0 x 10^^ that both of the methods 
are still poor in reproducing flattening, where the error bars for the results of the 
Fourier method are not shown in the figure, because they are too large. In the case 
of 6 = 2.0 X 10~^, on the other hand, a sharp bend at 9{ = Tr/2 and flattening at 
9 > ir/2 are clearly seen. Both of the methods provide consistent results within the 
errors. Here P{Q) is truncated at Q = 19 and only positive P{Q) are used. 

So far we have presented the results with only positive P{Q) in the singular f{6) 
model. Let us comment about the negative P{Q) in the model. For V = 150 in 
Fig. IHl for example, we employed positive P{Q) with Q = — 19. For larger values 
of Q, P{Q) oscillates and damps with a maximal amplitude of about 5 x 10~^. 
Negative values are also contained like P{Q) = 4.814 x 10"*, 2.910 x 10"*, -3.010 x 
10"*, -2.140 X 10"*, 2.696 x 10"*, 1.731 x 10"* and -2.395 x 10"* for Q = 18 - 24, 
respectively. By regarding these values as data and adding Gaussian noise to them, 
the resultant 2{9) in the Fourier method differs from that in Fig. |^ by at most 2 
% in the region 9{ < 9 < tt, while in the MEM about 6 % of maximal difference is 
found at 6* = 1.8(> 9{). 

§5. Summary 

We have considered lattice field theory with a 9 term, which suffers from the 
complex Boltzmann weight problem or the sign problem in numerical simulations. 
We have applied the MEM to the flattening phenomenon of the free energy originat- 
ing from the sign problem. In our previous analysis, we studied fictitious flattening of 
f{9). In the present paper, we have investigated whether the MEM could reproduce 
true flattening. For this, we have used mock data based on the simple models. In the 
flattening model, the MEM reproduces true flattening. This is different from what 
was obtained in (I)^^^ in the fictitious flattening case, where the MEM could yield 
monotonically decreasing behaviors of 2(9) in the whole region of ^(0 < 9 < vr). 
Whether the MEM could reproduce true flattening depends on the magnitude of the 
error. We have then investigated how the final images are affected by the errors in 
P{Q) by varying the value of 6. In the flattening model with 5 = 1/400, when pQ 
is as small as 1.0 x 10"'^, both of the results of the MEM and the Fourier method 
fail to reproduce flattening. However, as S decreases from 1/400, an immediate im- 
provement has been seen for the MEM. For a given value of po, the MEM becomes 
effective in reproducing correct images, when the value of 6 is reduced to the one 
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corresponding to the precision that X]q \^P{Q)\ ^ Po- With this precision of data, 
-2four(^) still receives much larger errors than those of the MEM. 

Since the flattening model is free from the singularity, we have considered a 
model, the singular f(6) model, which mimics a first order phase transition at 9{< vr). 
Although this is not a physical model in the strict sense, it would be useful for 
studying whether the MEM would reproduce the singular behavior. Such behaviors 
of the free energy at 9f are approximately reproduced depending on the precision of 
the data of P{Q). Here also, the MEM is effective when Y^q \5P{Q)\ < Z^. 

The errors of the final image in the MEM are calculated from the uncertainty 
of the image. In the procedure 3 in section 12.21 fluctuations of the image and the 
correlations among different 9 are properly taken into account according to the prob- 
ability. However, the smallness of the uncertainty in the MEM may not reflect correct 
magnitude of the errors. Estimate of such systematic errors is a task to be done in 
the future. On the other hand, the Fourier method reproduces the results reflecting 
the errors of P{Q) due to the propagation of the errors. This gives \5P{Q)\ as 
an estimate of the magnitude of the errors for any value of 9, which agree 

approximately with the results given in sections |21 and [l] In the limit of 5P{Q) 0, 
therefore, the Fourier method yields correct behavior, whose precision is only affected 
by the one of computations employed in the analysis; single, double or quadruple 
precision. 

In the present paper, we have concentrated on the issue whether the MEM would 
be effective with true flattening by employing mock data. A MEM analysis based 
on real Monte Carlo data will also be necessary. The results of such analysis for the 
QpAf-i j^Q^jg] ^[11 i^g presented in a forthcoming paper. 
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